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Abstract 

With the aim of finding evidence of tidal stripping of globular clusters (GCs) 
we analysed a sample of 13 elliptical galaxies taken from the ACS Virgo Cluster 
Survey (VCS). These galaxies belong to the main concentration of the Virgo 
cluster (VC) and present absolute magnitudes —18.5 < M z < —22.5. We used 
the public GC catalog of Jordan et al. (2008) and separated the GC population 
into metal poor (blue) and metal rich ( red) accord i ng to their integrated colors. 



The galaxy properties were taken from iPeng et al.l (120081 ). We found that 



1) The specific frequencies (SV) of total and blue GC populations increase 
as a function of the projected galaxy distances r p to M87. A similar result is 
observed when 3-dimensional distances r^o are used. The same behaviours are 
found if the analysis are made using the number of GCs per lO 9 Af (T). No 
correlations between Sn or T and r p or r 3 £, is observed for the red GC population. 
The correlations for the blue GCs (typically more extended) and the lack of 
correlations for the red GCs (more concentrated) with the clustocentric distance 
of the host galaxy are interpreted as evidence of GCs stripping due to tidal forces. 
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2) No correlation is found between the slope of GC density profiles of host 
galaxies and the galaxy distance to M87 (Virgo central galaxy). The lack of such 
a correlation is interpreted in terms of a shrinkage of the GC distribution after 
the stripping of GCs in the outermost region of galaxies. 

3) We also computed the local density of GCs (p ou t) located further than 
6.2 kpc from the galaxy center for nine galaxies of our sample. We find that the 
GC population around most of these galaxies is mainly composed of blue GCs. 
The two highest values of p out are found in the core of the VC (up to 100 kpc) 
and correspond to the two lowest values of Sn- 

Our results suggest that the number and the fraction of blue and red GCs 
observed in elliptical galaxies located near the centers of massive clusters, could 
be significantly different from the underlying GC population. These differences 
could be explained by tidal stripping effects that occur as galaxies approach the 
centers of clusters. 

Subject headings: galaxies: clusters: general — galaxies: elliptical and lenticular, 
cD — galaxies: interactions — Galaxy: globular clusters: general 



1. Introduction 

It has been demonstrated that tidal interactions can strongly affect the evolution of 
galaxies in clusters. Rapid gravitational encounters (galaxy harassment) as well as the global 
gravitational field of the cluster itself can dramatically change galaxy properties. Some of 
the properties that can be affected by tidal forces are those r elated to the population of 



globu lar clusters (GCs). Muzzio and collaborators in the 80s ( IMuzziol 11986 



Muzzio et al. 



19871 ). and more recently Bekki et al. (2003, hereafter B2003), used collisionless dynamical 
simulations to demonstrate that galaxies in clusters can lose a significant fraction of their GC 
population as a consequence of tidal stripping. Moreover, since the outermost regions of a 
galaxy are more sensitive to tidal effects, B2003 found that the number density profile of the 
GC population becomes steeper after the stripping. In addition, given that tidal encounters 
as well as the global effect of the cluster potential increase towards the inner region of the 
cluster, a positive correlation between the slope of the density profile of the GC population 
and the clustocentric distance of the galaxies is expected to be observed. 

If GCs are stripped from galaxies, a fraction of these GC should populate the intracluster 
medium. Bekki & Yahagi (2006, hereafter BK2006) investigated the spatial distribution 
of intracluster GCs (ICGC) using high-resolution cosmological simulations, and found that 
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ICGC population contribute 20-40% to the total GC population in high mass clusters, making 
the case for an inhomogeneous, asymmetric and elongated distribution. 



It is well known that color distributions of GCs are typically bimodal (jPeng et al.ll2006l ). 
indicating two sub-populations of GCs. Spectroscopic studies have shown that color bimodal- 
ity is mainly due to differences in the metallicity of the two sub-populations: the metal-poor 
GCs are blue, while the metal-rich GC are red. Blue GCs present a shallower and more 
extended radial density profile than red GCs, and consequently have a higher probability 
to be stripped from the parent galaxy. In the tidal stripping scen ario, the cluster co re and, 
therefore, the central galaxy play an important role. For example, iForte et al.l (119821 ) found 
that the properties of the GC population in M8 7 ("Virgo central ga laxy) are consistent with 
the accretion of GCs from less massive galaxies. ICote et al.l (119981 ) explored the possibility 
that the metal-poor population of giant elliptical galaxies arises from the capture of GC from 
other galaxies. Analysing the two brightest galaxies in Virgo, these authors concluded that 
it is possible to explain the observed properties of the blue GC population in M49 and M87 
without invoking new GC formation through mergers or multiple bursts. 



Forbes et al.l (119971 ) found evidence of tidal stripping in Fornax. These authors analyzed 



the specific frequency of GCs in four galaxies located in central cluster regions and found 
a marginal dependence of Sn (the number of GCs scaled to the galaxy luminosity) on the 
clustocentric galaxy distance r c . Nevertheless, the r c range consid e red is relatively small 
(r c < 200kpc). The same cluster was also analysed by iHilker et al.l (119991 ) . These authors 
suggested that an important fraction of the GC population of the cD galaxy NGC 1399 could 
be explained by GCs stripped from the central giant galaxies making those central galaxies 
to present low specific frecuencies. In order to quantify the importance of the tidal stripping 
effect, it is fundamental to study new clusters of galaxies and also to extend the limited 
range of clustocentric distances currently available. Recently, Peng et al. (2008, hereafter 
P2008) studied the formation efficiencies of GCs in the Virgo Cluster Catalog (VCC) sample. 
They found that both, high and low luminosity early type galaxies, have large Sn values, 
while those galaxies with intermediate luminosities (—20 < Mb < —17) present small and 
relatively constant Sn- They also found that Sn depends on the environment, i.e. 1) galaxies 
within a projected radius of 1 Mpc from the cluster center tend to have higher GC fractions, 
and 2) galaxies within r p ~ 40 kpc of M87 (and M49) have few or no GCs. The first result 
is interpreted in terms of the formation time, where GCs formation in central dwarf galaxies 
is biased because their stars form earlier; whereas the second result considers that GCs are 
tidally stripped by their giant neighbors. 

There are several processes that can produce tidal stripping of GCs. The main candi- 
dates are the pure gravitational forces due to the cluster potential well (see B2003), and the 
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galaxy -galaxy interaction (galaxy harassment). Using numerical simulations, iMoore et al. 
( 119961 ) demonstrated that in high mass galaxy cluster, this effect can indeed produce the 
stripping of an important part of t he stellar component o f galaxies. Based on the analysis 
of 228 elliptical galaxies in clusters, ICypriano et al.l (120061 ) found that galaxies in the inner 
region of clusters are 5% smaller than those in the outer regions. This result is interpreted 
as evidence of the stripping of stars from elliptical galaxies in the central regions of clusters. 

A precise knowledge of the number of GCs that galaxies could lose due to tidal stripping 
is particularly interesting for several reasons: i) If the models of GCs formation are going to 
be tested using observational data, GC population should be corrected by losses due to tidal 
stripping, ii) Some models for the formation and evolution of cD galaxies assume that a 



(Forte et al. 


1982. 


Hilker et al. 


1999) 



process will support this theory, iii) The GC population extends to large radii from the 
center of a galaxy. Therefore, GCs are ideal to study the outskirts of galaxies. If we can 
observe that GCs are being stripped by tidal interactions, this effect can be used to test a 
similar effect in other components of a galaxy, such as the dark matter halo. 

In this work, we study the GC population of a sample of 13 member galaxies of the 
Virgo Cluster (VC). The ACS Virgo Cluster Survey (ACS VCS) is described in Cote et 
al. (2004) (hereafter C2004). There are many rel evant paper s prod uced by the ACS team: 
Color distributi ons of the GCs are presented in iPeng et al.l (120061 ): GC size distributions 
are analysed in IJordan et al.l (120051 ): 3-dimensional distances for 84 galaxies in the VC are 
presented in lMei et al.l ( 120071 ): su rface brightness profile s, total magnitudes and colors of the 
sample galaxies are described in iFerrarese et al.l (120061 ): specific frequency, GC number, T 
parameter, and other galaxy parameters are presented i n P20 08, while the GC catalogue 
used in the present work is described in IJordan et al.l (120081 ) . The aim of this paper is 
to detect whether there exist observational evidence of the tidal stripping effect on the GC 
population due to galaxy-galaxy interactions and/or by the cluster potential well. The paper 
is organized as follows: sample selection is described in Section [2], we discuss the GC selection 
in Section 12.21 while in Section [31 we derive and analyze our results; the final conclusions 
are summarized in Section HI 
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2. The Sample 



2.1. Galaxy Sample 



The Virgo Cluste r (VC) is a well studied cluster at a distance of 1 6.5Mpc (ITonry et al. 



200lLlMei etaDl2007h . It has a velocity dispersio n of 776 ± 21 km s" 1 JGirardi et al.lll993h 



and an intracluster gas temperature of 2.4 keV (IDavid et al.lll993l ). Therefore, this cluster 
is an excellent candidate to study tidal stripping effects on the GC population. The ACS 
Virgo Cluster Survey is a dedicated optical and IR imaging survey of 100 Virgo's early- 
type galaxies, which made use of the Advanced Camera for Surveys (ACS) on-board the 
Hubble Space Telescope (HST). These observations provide a unique opportunity to study 
GCs in cluster galaxies, as done by the ACS VCS team. They also provide an insight into 
the GC populatio n in the intracluster medium (I CM). The complete survey is described in 



Cote et al.l (120041 ) and data reduction techniques in I Jordan et al.l (120041 ) . Briefly, the galaxies 



were observed through the F475W and F850LP filters that correspond to the Sloan g and z 
filters, respectively. The sample includes galaxies with B? < 16 and they were required to 
have morphological types of E, SO, dE, dE,N, dSO or dS0,N. 

The number of GCs depends on the luminosity of the host galaxy. Moreover, the 



specific frequency Sn depends on the galaxy magnitude (ILotz et al.ll2004l . iMiller fc Lotzll2007l . 
P2008). Since the ACS VCS sample spans over a wide luminosity range, and we want to avoid 
any luminosity bias in our analysis, we selected those galaxies with —22.5 < M z < —18.5. 
Over this luminosity range, the specific frequency Sn and the T parameter (defined in Sect. 
|3]) do not depend on the galaxy magnitude (P2008). On the other hand, since the GC 
population could depend on the host galaxy morphology, we have selected only elliptical 
galaxies. This type of galaxies shows (on average) the highest GC specific frequency for a 
given luminosity. In addition, we have only considered those elliptical galaxies that belong 
to the main cluster structure around M87 (IMei et al.l 120071 ). The last criterion excludes 
VCC1178 and VCC1025 since they belong to the cluster sub-structure associated to M49. 
The sample comprises 13 elliptical galaxies located around M87. They span a range of galaxy 
projected distance to M87 between 0.04 and 1.65 Mpc, while their total absolute magnitudes 
satisfy that —22.5 < M z < —18.5. Several galaxy properties, as the specific frecuency Sn, T 
parameter and other prop erties, were taken from P2008. We have also used the 3-dimensional 
distances r^o obtained by lMei et al.l (120071 ). For our sample we obtain clustercentric distances 
0.23 < r-i£, < 2.0. It shoul d be taken i nto a ccount that these distances have a line-of-sight 
uncertainties of ~0.6 Mpc. IMei et al.l (120071 ) adopted the standard ACDM cosmology with 
Q m h 2 = 0.127lQ;Q ( j ) 3 and h = 0.73 ± 0.03. Details of the 13 selected galaxies are quoted in 
Table [U 
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2.2. Globular Cluster Selection 



We specifically used the catalog of GC published by IJordan et al.l (120081 ). This catalog 
contains 12763 bona-fide GC candidates which have the probability pec > 0.5 of being 
a GC according to their selection procedure. Basically, the select ion procedure involves 
three main steps. The first one is described in IJordan et al.l (120041 ) and contemplates the 
modeling and removal of the galaxy s urface brightness distrib ution and subsequent object 
detection performed with SExtractor (IBertin fc Arnoutslll996l ). In the second step all the 
GC candida tes are run through the KINGPHOT code that fits PSF-convolved King models 
(King 1966) t o the ir surface brightness profiles. Details about this code can be found in 
Jordan et al.l (120051 ) . Evaluated parameters for each GC candidate are: magnitude z (total 
and in a 0.2arcsec diaphragm), celestial coordinates a and 5, half-light radius (rn) and 
concentration defined as the logarithm of the core and tidal radii ratio. However, after 
the initial selection of GC candidates mentioned above there are still residual contaminants 
such as foreground stars and background galaxies. These residuals are removed in the third 
step. At this stage, they consider the data as a mixture of points drawn mainly from two 
populations, namely the GC and the contaminants. Thus, the data can be modeled using 
a mixture model with two component, in which the total observed distribution in the r^-z 
plane is the result o f summing these two com ponents weighted by their respective sizes. The 
method follows thelFraley fc Ratfteryl (120021 ) procedure and it is also described in detail in 
Jordan et all (b()08h . 



3. Radial Dependencies 



3.1. Density Profiles of Globular Clusters 

The study of GCs around giant galaxies within clusters of galaxies has proved to be an 
useful tool for impr oving our knowledge about the formation and evolution of galaxies and 
their environments (IHanes et al.ll200ll ; ICote et al.ll200ll . 120031 ) . 



The radial distributi ons of GCs a r e often fitted wit 
as log(p) = alog(r) + (3 (IHarria Il986l . lAshman fc Zepi 



i a p ower law, which can be written 
1998 ). Typically, power law indices 



range from ~ —2 to ~ —2.5 for low-luminosity Es ([Harris & Harris! 120011 ). and are ~ —1.5 
for the most massive gEs. In general terms , blue metal-poor GCs present a shallower density 
profile than red GCs. iBassino et al.l (120061 ) found that both density profiles, for blue and red 
GCs, can be well fitted using log(p) = ar 1//4 + b. 

We have fitted a power law as well as r 1 / 4 law to our galaxy sample. GCs with g — z > 
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1.15 were considered metal rich (red GC), while GCs with g — z < 1.15 were considered 
metal poor (blue GC). The density profiles were computed using equal size bins (13.5 arcsec), 
therefore, the number of bins adopted in each galaxy depends on how far from the galaxy 
center we detected the GCs. In some cases only 3 bins were used, while in others up to 10 
bins were considered. Error bars were estimated using the bootstrap resampling technique. 
In Table [2] the fitting parameters for all the galaxies are quoted. The mean value obtained to 
the slope a for all GCs is a = —2. 13 ±0.08, while for red and blue GCs are a re d = —2.4 ±0.2 
and auue = — 1.7 ± 0.1, respectively. Particularly, Fig. [TJ shows the projected GC density 
profiles (r 1 / 4 law) for 4 galaxies of our sample. Left panels show the total GC density profiles, 
while right panels show the density profiles for the blue (squares) and red (triangles) GC 
populations separately. VCC1297 shows a flat density profile, therefore no profile was fitted, 
while VCC1422 does not present a red GC population. For the whole sample, we obtain a 
mean GC density profile with a = —1.59 ± 0.06, while for the red and blue GC populations 
the mean values are a re( i = — 1.8±0.1, and au U e = — 1.22±0.07, respectively. Seve ral galaxies 



in ou r sample do not present bimodality in the color distribution of GCs (see iPeng et al. 



20061 ). Nevertheless, for most galaxies in our sample, red GC density profiles are steeper and 
more centrally concentrated than the corresponding blue GC density profiles. This could 
indicate that we indeed have two different GC populations. 

Towards the inner region of the galaxy cluster, tidal effects are expected to become 
stronger due to a higher probability of encounters among galaxies, as well as the increase of 
the potential well. If we consider, in addition, that GCs located in the outskirts of galaxies 
are more prone to be stripped, then a correlation between the slope of the GC density profiles 
and the distance of hosts galaxies to M87 should be observed. This effect has already been 
reported by B2003 using numerical simulations. Nevertheless, for both, blue and red GCs, 
we do not find any correlation between these two parameters (see Fig. [2]). A possible 
interpretation to this result is given in Sec. HI 



3.2. Specific Frequency Sn and T parameter 



Harris fc van den Berghl (119811 ) introduced the specific frequency Sn as a measure of 

(1) 



GC richness normalized to the host galaxy luminosity: 

S N = N GC 10 0A ^ +1 ^ 



where Nqc is the total number of GCs and M y is the V^-ban d absolute magnitude. Typical 
values of Sn are ~1 for spirals a nd SO galaxies (jBarmbyi 120031 ). ~3 for ellipticals, and ~6 for 
cD galaxies (jOstrov et al.lll998l ). Nevertheless, the correlation between morphological types 
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and Sn contains some dispersion. On the other hand, several early studies suggest that Sn 
increases with galaxy luminosity (iMiller et al.lll998l . lLotz et al.l l200ll . iMiller fc Lotzl 120071 ) . 
However, P2008 found that the GC mass fraction is high i n both giants and dw arfs galaxies, 
being universal in intermediate luminosity galaxies, while IStrader et al.l (120061 ) only found a 
weak signal in their Sn,b — Mb data. 



Sn can also depend on the environment, iFerguson fc Sandagd (Il989l ) found that the 
number density of nucleated dwarf elliptical galaxies (dE) is more centrally con centrated than 
that o f bright non-nucleated dE galaxies in the Virgo and Fornax Clusters. IMiller fc Lotz 
(120071 ) did not find a clear dependence of Sn on the environment for dEs. Using the VCS, 
P2008 found that galaxies within a projected radius of IMpc from M87 tend to have higher 
GC fractions, and this effect was interpreted in terms of the galaxy formation time. It should 
be noted that the described effect contradicts the tidal stripping scenario. Nevertheless, this 
result is observed when using the total sample of galaxies and does not discriminate by their 
luminosity, morphology or location Virgo's sub-structures. 

Panels (a) and (b) of Fig. [3J show Sn as a function of (r p ) and (r 3D ) respectively 
for our sample of elliptical galaxies (see Table [T]). We observe that larger Sn values are 
found at larger galaxy distances to M87 (~ Virgo center). This result is indicating that 
the normalized number of GCs decreases as the host galaxy is closer to the Virgo center, 
where the tidal stripping is expected to be more efficient. A similar behavior is observed 
in panels (c) and (d), where we plot the specific frequency in the z bandpass Sn,z, as a 
function of the galaxy distances to M87. In the aforementioned Figure, panels (e),(f) and 
(g),(h) show the specific frequency SN,z,Uue and SN, z ,red versus galaxy distances to M87 for 
the blue and the red GC populations, respectively. It can be seen that both correlations, 
projected Sn vs. r p and Sn vs. r 3D , are only present for the blue GC population. The 
lack of correlation for the red GCs could be explained by the fact that the distribution of 
red GC is more highly concentrated than the equivalent blue GC distribution, and therefore 
less affected by tidal effects. In order to quantify this effect, we split our sample into two: 
galaxies with r p < 0.8Mpc and r p > 0.8Mpc. The adopted threshold for r 3 r> is 1.3Mpc 
(assuming that r 3 D = y/3r p ). For each sub-sample, the mean Sn are shown in Table El As it 
can be seen, galaxies in the outskirts of the Virgo cluster have mean Sn typically larger than 
those calculated for galaxies closer to M87. Particularly, we excluded VCC 1475 to compute 
the mean SN, z ,biue and the mean Tu U e due to its unexpected high values. Alternatively, for 
the blue GCs we fitted a linear regression between SN, z ,Uue and r p and r 3 o (see panels (e) 
and (f) in Fig. [3]). Table @] shows the fitting parameters, where it can appreciated that fitted 
slopes differ from zero at six sigma levels. For the fits, we also excluded VCC 1475. 



In addition, we show in Fig. [3j Sn vs. r p (left panels) and Sn vs. r 3 D (right panels), 
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discriminating the galaxy sample according to the absolute magnitude M z of the galaxies. 
It can be appreciated that the correlation between Sn and r p and r 3D is independent of the 
galaxy luminosity. 

Usin g stellar masses o f gala xies in the VCS, P2008 calculated the T parameter intro- 
duced by lZepf fc Ashman! (119931 ). defined as the number of GCs per 1O 9 A^0, 

T = — (2) 

where A4g is the stellar mass of the host gal axy. The use o f T all ows the comparison across 
galaxies with different mass-to-light-radius. iMiller fc Lotzl (120071 ) and P2008 found that T 
correlates with galaxy luminosity, in the sense that higher T values correspond to brighter 
galaxies. Nevertheless, in the magnitude range of our sample, no correlation between T and 
the absolute magnitude is observed. Panels (a) and (b) of Fig. H] show T vs. r p and T vs. 
r 3 £i, respectively. Panels (c) and (d) correspond to the blue GC population for which we 
fitted a linear regression (again, the r 3D fit excludes VCC 1475), whereas panels (e) and (f) 
show the corresponding correlations for the red GCs. As it was observed for Sn, we only 
find correlations for the blue and the total GC populations. Higher values of T are found in 
the outskirts of the Virgo cluster (see also Tab. [3] and Tab. @J. 



3.3. Local background of GCs 

The observed population of GCs is superposed on a background contamination that 
exponentially grows with the apparent magnitude limit. This contamination is basically 
constituted by high redshift galaxies, foreground stars and a projected ICGC population. 
If tidal stripping of GCs is an efficient mechanism in clusters of galaxies, a fraction of the 
ICGC background would also be formed by stripped GCs. Moreover, this ICGC background 
would show a gradient, having higher values towards smaller projected clustocentric dis- 



tances. There are several works that suggest the existence of the ICGC. IWest et al.l (119951 ) 
suggested that the variation in the GC population in supergiant elliptical galaxies in clusters 
points towards the existence of a GC population that is not bound to a ny galaxy and moves 



throughout the core of the galaxy cluster. This was also suggested by iJordan et al.l (120031 ) 
who used HST images of the rich cluster A1185. This galaxy cluster contains no bright 
galaxies, and shows an excess of point sources that is consistent with the existence of an 
ICGC population. The efficiency of GC tidal stripping due to the cluster potential well can 
be evaluated in terms of the tidal radii. Given a galaxy with mass A4 ga i and at a given 
distance r c from the cluster center, the tidal radius can be computed as 

rtuu = [GM gol /2vl rc ] x '^ (3) 
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where v C i rc is the circular velocity of the galaxy cluster potential. Th erefore, for a typical 
elliptical galaxy with a total mass of 10 12 A4 & (jSchindler et al.lll999l ). located in the VC 



potential, we obtain tidal radii of 11, 14 and 17 kpc for clustocentric distances of 0.2, 0.5 
and 1.0 Mpc, respectively. These numbers indicate that if a galaxy passes close to the VC 
center it could lose part of its GC population due to tidal stripping. 

Numerical simulations run by B2003 and BK2006 suggest that ICGCs are not necessarily 
uniformly distributed. GCs stripped from galaxies have not had enough time to be uniformly 
distributed in the IC medium. In fact, the simulation of B2003 shows that a fraction of 
the stripped GCs produces an enhancement in the local density around the parent galaxy. 
Moreover, BK2006 showed that, on average, ICGCs appear to have projected radial density 
profiles that grow towards the galaxy cluster centers. Based on these two results, we study the 
possible dependence of the ICGC population on the galaxy distance to M87. In our sample, 
we find that galaxies that are fainter than M z = —20.8 (9 out of 13 galaxies) have a flat 
GC projected density profile for galactocentric distances greater than 1.3arcmin (6.2 kpc). 
These galaxies are suitable to determine the local GC density around them since they have 
relatively small apparent diameters compared to the ACS field. We then calculated the local 
GC density beyond 6.2 kpc from the center of each of these galaxies. However, this density 
is still contaminated with high redshift galaxies and foreground stars. Hereafter we will refer 
to this local density as p out . 

P2008 determined fore- and background contamination using 10 blank fields images 
taken with the Wide Field Planetary Camera 2 (WFPC2). They found that there are on 
average 5.8 ± 2.5 contaminants sources per WFPC2 field. We use this number to estimate 
the contribution of the mean fore- and background contamination to p out . Figure [5] shows the 
total, blue and red (from top to bottom) local density p out around the selected nine galaxies 
as function of r p (left panels) and r 3 £> (right panels). Error bars were estimated using the 
bootstrap resampling technique. Horizontal dashed lines correspond to the mean value of 
p ou t computed for those galaxies with clustocentric distances r p > 100 kpc. These values are 
2.4±0.3 (arcmin- 2 ) for all GCs, 1.7±0.3 for blue GCs, and 0.6±0.2 for red GCs. Horizontal 
dotted lines show the mean density (blue and red) of the fore- and background contaminants 
derived by P2008 from 10 blank fields (1.1 ± 0.5). It can be seen that nine (eight) galaxies 
have values of p out for the total (blue) GC population consistent or higher than expected 
mean values of the fore- and background contaminants. Except for two galaxies, the values 
of p ou t for the red GC population, are lower than the expected for the total mean fore- 
and background contamination. On the other hand, we found for the total and blue GC 
populations, that three galaxies (VCC1297, VCC1279 and VCC0828) have p ou t values three 
times larger than the contaminants mean density. The same figure also shows that the two 
highest values of p ou t correspond to the galaxies with the lowest values of Sn (VCC1297 and 
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VCC1279) that are also the galaxies closest to M87 (see Tabled]), which has an extended halo 
of GCs. The extended GC halo of M87 could be partially contributing to the high density 
levels observed around these two galaxies. However, since M87 is located near the center of 
the cluster potential well we are unable to determine if this additional contribution is derived 
from GCs bound to M87 itself, or if they are truly ICGCs. In any case, the procedure we use 
to compute the background is independent of the origin of GCs and therefore, our method 
does accurately represent the GC background around these galaxies. 



4. Conclusions and Discussion 

We analysed a sample of 13 elliptical galaxies taken from the ACS VCS with the aim 
of finding evidence of the tidal stripping of GCs. Particularly, we focused our attention 
on the dependence of several galaxy properties with the radial clustocentric distance. The 
selected sample comprises host galaxies in the magnitude range —22.5 < M z < —18.5 that 
belo ng to the main conc entration of the VC. The GCs around these galaxies were taken 



from I Jordan et al.l (120081 ) . while several galaxy properties were taken from P2008. The color 



index g — z was used to discriminate metal poor (blue) from metal rich (red) GCs. 

There is a correlation between the specific frequency, in the V and z-bands, and the 
galaxy distances to M87. Larger values of Sn are found at larger clustocentric distances. 
The effect is observed using both projected or 3-dimensional distances. If GCs are selected 
according to their color, the correlation is only observed in the blue (metal poor) population. 
These results are in agreement with the idea that those galaxies closer to the Virgo center 
are loosing a fraction of their GCs due to tidal stripping. The lack of correlation for the red 
GC (metal rich) population could be explained if we consider that these objects have a more 
concentrated distribution around the host galaxy, and therefore they are less affected by 
tidal effects than the blue GCs. Similarly, we found a correlation between the T parameter 
(the number of GCs normalized to galaxy stellar mass) and r p and T and r 3 £>, for the total 
and blue GC populations only. T increases as the galaxy distance to M87 increases. The 
lowest values of Sn and T are found within lOOkpc to the cluster center, where P2008 have 
also reported evidences of tidal stripping. 

As it was explained before, P2008 found that: a) Sn and T depend on the environment, 
i.e. galaxies within a projected radius of 1 Mpc from the cluster center tend to have higher 
GC fractions than galaxies at larger clustocentric distances, b) Galaxies within r p ~ 40 kpc 
from M87 (and M49) have a few or no GCs. The first result is interpreted in terms of 
formation time. GC formation in centrally located low-mass galaxies is biased because their 
stars form earlier. The second result considers that GCs are tidally stripped from the parent 
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galaxies by their giant neighbors. At first glance, our findings appear to differ from the 
aforementioned P2008 result (a), since we find that the mean Sn value is indeed higher 
at larger r p . Nevertheless, before we can compare our results to those of P2008 we have 
to consider the following: i) although all VCC galaxies are early-types, P2008 analysed the 
radial dependence for a great variety of galaxy morphologies, while our analyses only consider 
pure elliptical galaxies; ii) our sample excludes bright and faint VCC galaxies and comprises 
only those objects in the magnitude range —22.5 < M z < —18.5, for which P2008 data show 
no correlation between the galaxy magnitude and S^. This is the most important difference 
between our analysis and that performed by P2008, since the dependences of Sn and T 
on the galaxy magnitude tend to erase any other correlation; iii) our study only considers 
galaxies that belong to the main concentration around M87 and excludes galaxies close to 
M49. It is more likely that the dynamic history of these galaxies is more related to the group 
of M49 than to the Virgo cluster. 

We tested a possible correlation between the slopes of the GC density profiles and the 
clustocentric distances as predicted by B2003 numerical simulations . Computed slopes for 



our galaxy sample do not show correlation with r p nor with r 3D . ICypriano et al.l (120061 ) 
found that galaxies in the inner regions of clusters are 5% smaller than t hose in the outer 
regions. Th ese authors interpret this result in terms of star stripping. lAguilar &: White 



( 19861 . Il987l ) studied this effect using numerical simulations and they found that the r 1//4 
luminosity profile is robust, since it is recovered after the galaxy interaction or merger. 
However, r 1//4 paramenters, the effective radius r e and the surface brightness p e inside r e , 
change depending on the collision type. Strong collisions result in a reduced r e and a higher 
/i e , while weak collisions produce the opposite effect. If GCs are affected by the same process 
that affect the stars, the GC distribution could also be modified, but the final density profile 
will also depend on the interaction class. If this model is correct, no correlation between the 
slopes of the GC distribution and the clustocentric distances should be expected. 

We also computed the GC density p out , outside 6.2 kpc from the galaxy center for nine 
galaxies of our sample. Our results show that values of p out for the total (blue) GC population 
around nine (eight) of the selected fields are consistent, or larger, than those expected from 
the fore and background contaminants. Except for two of the nine selected fields, the values 
of Pout for the red GC population are lower than expected for the mean fore- and background 
contamination. These results indicate that the ICGC population around most of the galaxies 
of our sample is mainly composed of blue GC. The two highest values of p ou t (that correspond 
to the two lowest values of Sn) are found in the core of the VC (within lOOkpc), confirming 
P2008 results. They could be stripped GCs from neighbor galaxies, some fraction of the 
very extended GC distribution of M87, and/or a genuine ICGC population. The two fields 
at larger clustocentric distances with total and blue values of p ou t larger than expected for 
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the fore- and background contamination could be either a cosmic variance or the result of 
GCs stripped from the parent galaxy that have not had enough time to uniformly distribute 
into the intracluster medium. 

The evidence for tidal stripping found in this work suggests that special care should be 
taken in studies that attempt to model the formation and evolution of the GC population. 
The number, proportion and distribution of GCs observed in those elliptical galaxies close 
to the center in a massive cluster, appear to significantly differ from the GC population 
in galaxies before they approached the center of the cluster. Moreover, a good estimate of 
the background contamination in GC counts is critical, since we have shown the possible 
existence of both, a mean gradient as a function of r p , and a possible cosmic scatter of 
the mean local background around some galaxies. Finally, more realistic simulations are 
needed in order to confront observational results with model predictions. A complete set of 
cluster galaxies selected from a cosmological simulation should be analyzed. Besides, metal 
poor and metal rich GCs with their corresponding radial distributions should be considered. 
Simulations constrained in this way, have the potential to predict the initial undisturbed 
radial distribution of GCs around galaxies in rich clusters. 
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Fig. 1. — GC projected density profiles for four galaxies in our sample. Filled squares, 
open triangles and open squares correspond to the total, red and blue GC, respectively. 
Long dashed lines and short dashed lines correspond to r 1 / 4 fit for the blue and the red 
GC population, respectively. Errors bars were estimated using the bootstrap resampling 
technique. 
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Fig. 2. — Slopes of the GC projected density profiles as a function of the galaxy projected 
distance r p (left panels) and as a funtion of the galaxy 3-dimensional distance r^r> to M87 
(right panels). From top to bottom: correlations for the total GC population, blue and red 
GC pupulations. 
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Fig. 3. — S N of each galaxy as a function of the clustocentric projected distance r p of the 
host galaxy (left panels) and as a funtion of the galaxy 3- dimensional distance r S D (right 
panels). From top to bottom: correlations for the total GC population in V-band, total GC 
population in the z-band, blue GC pupulation in the z-band and red GC population in the 
z-band. Black solid lines in panels e) and f) correspond to linear regressions. Open and solid 
symbols consider two different ranges in the absolute magnitude M z of the host galaxy. 
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Fig. 4. — Some as Fig. [2] but now for the T parameter. 
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Fig. 5. — Total local density p out as a function of the distance to M87, from top to bottom: 
for the total GC population, blue GC pupulation and red GC population. Black solid lines 
in panels e) and f) correspond to linear regressions. The horizontal dashed lines correspond 
to the mean total density computed around galaxies at r p > 100 kpc (vertical lines). Dotted 
lines show the expected mean density of the fore- and background contaminants. Errors bars 
were estimated using the bootstrap resampling technique. 



Table 1. Galaxy properties in ACS Virgo Cluster Survey taken from lPeng et al.l (120081 ) and lMei et al.l (120071 ). 



vcc 


M, 




r p 




>"3D 




N G C 












■5 N, 






S 










d 


T 












Type 


(!) 


(2) 




(3) 




(4) 






(5) 




(6) 








C) 








(8) 






(9) 






(10) 


(11) 






(12) 




(13) 


1903 


-22 


19 


0.82 


1 


87 ± 





35 


803 ± 355 


2 


71 ± 


1 


19 


1 


07 ± 





47 





55 ± 


38 





52 ± 





27 


9.6 ± 4.8 


6.2 ± 2 





4 


1 ± 1 


4 


B5 


1231 


-21 


58 


0.31 


1 


30 ± 





36 


376 ± 97 


1 


98 ± 





51 





88 ± 





23 





50 ± 


IS 





38 ± 





13 


7.0 ± 2.4 


4.0 ± 1 


7 


3 


± 1 


3 


B5 


2000 


-20 


66 


1.03 


2 


03 ± 





35 


205 ± 28 


2 


97 ± 





41 


1 


12 ± 





15 





97 ± 


13 





15 ± 





08 


8.6 ± 1.9 


7.5 ± 1 


6 


1 


2 ± 


7 


E5 


1664 


-20 


97 


0.48 





77 ± 





38 


213 ± 31 


2 


35 ± 





34 





87 ± 





13 





54 ± 


10 





33 ± 





07 


8.2 ± 2.3 


5.1 ± 1 


5 


3 


1 ± 1 





Eli 


1279 


-20 


72 


0.04 





47 ± 





40 


58 ± 11 





72 ± 





15 





30 ± 





06 





21 ± 


15 





09 ± 





04 


2.7 ± 0.7 


1.8 ± 


(i 





8 ± 


3 


E2 


1619 


-20 


19 


0.33 


1 


51 ± 





36 


84 ± 19 


1 


52 ± 





34 





70 ± 





l(i 





52 ± 


13 





19 ± 





09 


4.9 ± 1.4 


3.6 ± 1 


1 


1 


3 ± 


7 


E/SO 


0828 


-20 


07 


0.38 


1 


42 ± 





42 


69.5 ± 9.8 


1 


52 ± 





21 





65 ± 





09 





50 ± 


08 





50 ± 





08 


5.1 ± 1.3 


3.9 ± 1 


1 


1 


2 ± 


4 


E5 


1630 


-22 


19 


0.82 


1 


87 ± 





35 


47 ± 11 


1 


11 ± 





26 





44 ± 





10 





26 ± 


08 





18 ± 





06 


4.0 ± 1.3 


2.4 ± 


9 


1 


7 ± 


7 


E 


1146 


-20 


06 


0.34 





51 ± 





38 


72 ± 12 


1 


93 ± 





32 


1 


00 ± 





17 





30 ± 


12 





70 ± 





12 


8.3 ± 2.5 


2.5 ± 1 


1 


5 


8 ± 1 


8 


E 


1913 


-19 


74 


1.58 


1 


67 ± 





40 


71 ± 14 


2 


02 ± 





40 





90 ± 





IS 





74 ± 


15 





16 ± 





09 


6.6 ± 2.1 


5.4 ± 1 


8 


1 


2 ± 


8 


E 


1475 


-19 


42 


1.13 


1 


15 ± 





38 


81 ± 10 


3 


05 ± 





40 


1 


38 ± 





IS 


1 


33 ± 


17 





06 ± 





06 


10.5 ± 2.5 


10.1 ± 2.4 





4 ± 


5 


E 


1422 


-18 


73 


0.62 


1 


60 ± 





43 


24.9 ± 6.0 


1 


62 ± 





39 





80 ± 





10 





54 ± 


l(i 





26 ± 





11 


6.5 ± 2.8 


4.4 ± 2.0 


2 


1 ± 1 


2 


E 


1297 


-18 


75 


0.04 





23 ± 





46 


4 ± 11 





42 ± 


1 


00 





16 ± 





37 





11 ± 


30 





05 ± 





22 


1.2 ± 3.0 


0.9 ± 2.4 





4 ± 1 


7 


E 



Note. — (1) Number in Virgo Cluster Catalog. (2) Absolute z magnitude. (3) Projected distance from M87, in Mpc. (4) 3-dimcnsional distance from M87, in Mpc. 
(5) Total number of GCs. (6) Specific frequency in V-band. (7) Specific frccucncy in z bandpass. (8) Specific frccuency in z for blue GCs. (9) Specific frccucncy in z 
for red GCs. (10) T parameter: Nqq normalized to stellar mass of IO^Mq. (11) T parameter to blue GCs. (12) T parameter to red GCs. (13) Morphological type. 
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Table 2. Fit parameters of density profiles. 
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Table 3. Mean values of Sn 


and T as 


a function of 


r p and r 3D . 
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In the r$r> analysis, we excluded VCC1475 (sec the text). 
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Table 4. Parameters of the linear fit of SN, z ,biue and Tu U e as function of r p and r^D- a is 

the slope and b the zeropoint. 
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